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We report the results of a search for the rare baryonic decays B° — ► pp, A A, and B + — > pA. The 
analysis is based on a data set of 31.7 x 10 6 BB events collected by the Belle detector at the KEKB 
e + e~ collider. No statistically significant signals are found, and we set branching fraction upper 
limits B(B° -»■ pp) < 1.2 x 1(T 6 , B(B° -> A A) < 1.0 x 1(T 6 , and B(B+ -> pA) < 2.2 x 10~ 6 at the 
90% confidence level. 
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FIG. 1. Illustrative diagrams for B decays to charmless 
baryon pairs. 

The Belle collaboration recently reported the observa- 
tion of the decay process B + — ► ppK + |l) , which is the 
first known example of a B meson decay to a charmless 
final state containing baryons. In this paper we report 
the results of a search for the related two-body modes 
B" -> pp, AA and B+ -> pA §. In the Standard Model, 
these decays are expected to proceed via color-suppressed 
b — > u tree diagrams (Figs. 1(a) and (c)) and b —> s, d 
penguin diagrams (Figs. 1(b) and (d)). The search is 
based on a 29.4 fb _1 sample of e + e~ data accumulated 
at the T(4S I ) resonance, which contains 31.7 million BB 
pairs. A previous search for these decays by the CLEO 
collaboration using a 5.41 fb _1 sample of T(45) data 
yielded 90% confidence-level (C.L.) upper limits M. 

Belle 0j is a general purpose detector operating at 
the KEKB asymmetric e + e~ collider ||. Tracking in- 
formation is provided by a silicon vertex detector and 
a central drift chamber (CDC) in a 1.5 Tesla magnetic 
field. Hadron identification (PID) for ir/K/p discrimina- 
tion is obtained from CDC dE/dx measurements, aerogel 
Cerenkov counter pulse heights, and timing information 
from the time-of-flight system. Electron identification is 
based mainly on CsI(Tl) electromagnetic calorimeter and 
CDC dE/dx information. Kl and muons are identified 
by a system of resistive plate counters interleaved with 
the iron plates of the flux-return iron yoke. 

The event selection criteria are based on tracking and 
PID requirements, and are optimized using Monte Carlo 
(MC) simulated event samples. 

All primary charged tracks are required to satisfy the 
following track quality criteria based on the track impact 
parameters relative to the interaction point (IP), which 
is determined run-by-run. The z axis is defined by the 
positron beam line. The deviations from the IP position 
are required to be within ±0.05 cm in the transverse (x-y) 
plane, and within ±2 cm in the z direction. Tracks that 
satisfy the muon or electron identification requirements 
are rejected. 
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pjc" Mass (GeV/c 2 ) 
FIG. 2. The mass distribution of the selected A — > pn~ 
candidates for a typical run period. The MC distribution is 
shown as a dashed histogram. The ±5 MeV mass window is 
indicated by the vertical arrows. 

Primary proton candidates are selected based on 
p/K/ir likelihood functions obtained from the hadron 
identification system. The selection criteria are L _^ Lk > 

0.6 and L L ^ L > 0.6, where L p / K /^ stands for the pro- 
ton/kaon/pion likelihood. 

A candidates are reconstructed via the pir~ decay 
channel and are selected using cuts on four parameters: 
the angular difference between the A flight direction and 
the direction pointing from IP to the decay vertex in the 
transverse plane; the distance between each track and 
the IP in the transverse plane; the distance between the 
decay vertex and the IP in the transverse plane; and 
the displacement in z of the closest approach points of 
the two tracks to the beam axis. The secondary protons 
are required to have L L ^ L > 0.6. The pir~ mass spec- 
trum after the application of the above selection criteria 
is shown in Fig. 2 for a typical run period. The peak po- 
sition is consistent with the nominal A mass H and the 
mass resolution is about 0.9 MeV/c 2 . Finally, we require 
the invariant mass of the A candidate to be within ±5 
MeV/c 2 of the nominal A mass. 

Due to the high momentum of the primary particles 
from these two-body decay modes, the background from 
generic B decays is negligible. This is checked with MC 
samples of B + B~ and B°B° pairs where the B-mesons 
decay dominantly via b — > c processes. We also checked 
backgrounds using MC samples of charmless B meson 
decays that include low multiplicity B decays into final 
states with n, K, K*, p, u>, </>, ry, and rf mesons. Only 
the pp mode is found to have any background contamina- 
tion. However, the level of this contamination is negligi- 
ble, corresponding to less than one event for the current 
data set. 

The main background is from continuum qq processes. 
This is confirmed using an off-resonance data set (2.3 
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TABLE I. Results of search for the exclusive baryon modes. The signal yields, Y, and errors are determined from maximum 
likelihood fits. The 90% C.L. upper limits from the fits and from the counting method are listed together. We quote the higher 
values as our conservative estimates for upper limits. The efficiencies, e, are obtained from MC. The 90% C.L. upper limits for 
the branching fractions, B, determined by this experiment are shown along with previous CLEO results. 



Mode 


Y 


UL (fitting/counting) 


e(%) 


S(HT 6 ) 


CLEO S(10~ 6 ) 


PP 


6 +a7 

u -°-0.6 


7.0/9.7 


27.5 ± 2.0 


< 1.2 


< 7.0 


AA 


o o +0 - 5 

u - u -o.o 


3.0/3.2 


10.8 ± 1.1 


< 1.0 


< 3.9 


pA 


1 o+ 2 - 5 

L - u -i.o 


7.0/10.4 


16.2 ± 1.4 


< 2.2 


< 2.6 



ft) -1 taken 60 MeV below the T(45)) and a MC sam- 
ple of 65 million continuum events. These continuum 
events have a jet-like topology while BB events are more 
spherical in the T(45) center of mass (CM) frame. For 
continuum event rejection we use cos 9t, the cosine of the 
angle between the direction of one primary decay particle 
and the thrust axis |?| of the non-candidate tracks and 
showers. This distribution is nearly flat for signal events 
and is strongly peaked at ±1 for continuum background. 
We also use cos 6b , the cosine of the angle between the 
B candidate flight direction and the positron beam di- 
rection. The signal has a sin 2 6b distribution while the 
background is uniformly distributed. We require the ab- 
solute values of cos#t and cos6 B to be less than 0.9 for 
AA decays and less than 0.8 for the other modes. For 
the latter case, the background reduction factor is more 
than five, while ~70% of the signal is retained. 

We use the following two kinematic variables to iden- 
tify the reconstructed B meson candidates: the beam 
constrained mass, m bc — E^ cam — p B , and the energy 
difference, AE = E B - E hc&m , where Ehcam, Pb and E B 
are the beam energy, the momentum and energy of the 
reconstructed B meson in the T(4S*) CM frame, respec- 
tively. We retain events with 5.20 GeV/c 2 < mbc < 5.29 
GeV/c 2 and -0.2 GeV < AE < 0.2 GeV. The signal 
yield is extracted by maximizing the likelihood function 

N 

L = e-^+V l[[sP s {m hCi , AEi) + bP b (m hCi> AEi)}, 

where N is the total number of candidate events, s(b) de- 
notes the signal (background) yield and P s ^) denotes the 
signal (background) probability density function (pdf). 

The mbc and AE signal pdf's are determined by MC. 
We use a Gaussian function as the signal pdf for the 
TObc distribution and a sum of two Gaussians for the AE 
distribution. The Gaussian parameters (mean and a) 
are determined separately for each mode. Background 
shapes are determined from events in sideband regions of 
AE and mbc separately. We adopt an empirical function 
Q to model the mb c background shape (for events with 
0.1 GeV < \AE\ < 0.2 GeV) and a first-order polyno- 
mial for the AE background shape (for events with 5.20 
GeV/c 2 < m bc < 5.26 GeV/c 2 ). 



Table 1 summarizes the results. The efficiencies for the 
AA and pA modes include the A — > pn~ branching frac- 
tion (64%). The efficiencies are determined from signal 
MC with the identical event selection and fitting proce- 
dure as for the data. Figures 3 and 4 show the mbc (with 
| AS | < 0.05 GeV) and AE (with 5.27 GeV/c 2 < m bc < 
5.29 GeV/c 2 ) projections for these three modes, respec- 
tively. Projections of the fits are shown as smooth curves. 
No statistically significant signals are found and we de- 
termine 90% C.L. upper limits on the signal yields by in- 
tegrating the likelihood function. We also compute limits 
using a counting method We define a signal region 
by 5.27 GeV/c 2 < m bc < 5.29 GeV/c 2 and \AE\ < 0.05 
GeV, and treat the number of background events from 
the maximum likelihood fit as a prediction for the back- 
ground in this region. We then count the number of 
events actually observed, and apply the method of |jj to 
obtain 90% C.L. intervals. 

To estimate the possible influence of fluctuations on 
the determination of the upper limits, we vary the param- 
eters of the pdf's by one standard deviation and change 
the form of the AE signal pdf to a single Gaussian. The 
relative change of the upper limit is 25%, which is mainly 
due to the changes in the background shape. The upper 
limit determination by the counting method is checked by 
redefining the signal region (varying the AE range from 
2a to 4<r) and comparing the outcomes. The upper limits 
from both methods listed in column 3 of Table 1 include 
these fluctuations. We quote the values from the count- 
ing method as the most conservative upper limits. Fit 
projections with the signal yield fixed at the upper limit 
for each mode are shown as the superimposed dashed 
curves in Figs. 3 and 4 for comparison. 

The systematic error due to the efficiency of the pro- 
ton identification (p-K and p-ir) criteria is studied using 
A samples. We vary the likelihood ratio requirement for 
protons and compare the ratio of reconstructed A yields 
in data and MC. The overall error is about 3%. We in- 
clude a 2% error per track to account for the uncertainty 
in tracking efficiency. The A reconstruction efficiency is 
checked by comparing the flight distance distributions of 
data and MC. They agree very well and no additional 
error is assigned. The correlated parts of the errors are 
added linearly to obtain the overall uncertainty in the 
tracking efficiency and the uncertainty in the P1D effi- 
ciencies (for p and p), then the resulting errors are com- 
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m bc (GeV/c 2 ) 



FIG. 3. The distributions of m bc for (a) B° -> pp, (b) 
B° — > AA and (c) B + —* pA candidates. The solid curve 
is the projection of the maximum likelihood fit. The dashed 
curve shows the fit with the signal yield fixed at the 90% CL. 
upper limit. 
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FIG. 4. The distributions of AE for (a) B° -» pp, (b) 
B° — > AA and (c) _B + — » pA candidates. The solid curve 
is the projection of the maximum likelihood fit. The dashed 
curve shows the fit with the signal yield fixed at the 90% C.L. 
upper limit. 



bined in quadrature. When determining the upper limit 
for the branching fraction, the efficiency was reduced by 
one standard deviation. The efficiencies and upper limits 
for all three decay modes are listed in Table 1. 

In summary, we have performed a search for the rare 
baryonic decays B° — > pp, AA, and B + — > pA with 31.7 
million BB events collected by the Belle detector at the 
KEKB e + e~ collider. No statistically significant signals 
are found for these modes, and we set upper limits on 
their branching fractions at the 90% C.L.. The upper 
limits are: 



B(B° 
B(B° 
B{B+ 



>pp) < 1.2 x 1CT 6 ; 
AA) < 1.0 x 10~ 6 ; 
■pA) < 2.2 x 10~ 6 . 



These are currently the most stringent limits for these 
decays. The limit on B — > pp has been improved by a 
factor of six compared to the existing bound ||. These 
limits are considerably below corresponding charmless 
mesonic branching fractions which are typically at the 
level of 10 -5 . In contrast, the recent Belle measurement 
of B + — > ppK + jlj indicates a relatively large branch- 
ing fraction of 4.3 x 10~ 6 . The underlying physics ]Tc| ] 
is now being probed with rapidly accumulating data and 
new experimental results. 
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